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Within the framework of mode-coupling theory, the glass-transition scenario is inves- 
tigated for a system of particles interacting with a hard-core repulsion and an additional 
square-shoulder soft core at larger distances. The static structure is calculated from the 
potential in Percus-Yevick approximation. For certain widths of the shoulder, the exponent 
parameter A along the glass-transition line shows a double peak. At both peaks, A can reach 
unity indicating the existence of higher-order glass-transition singularities. It is shown that 
these higher-order singularities originate from a line of avoided glass-glass transitions. 

§1. Introduction 

In the field of glassy slow dynamics, many experiments and simulations have been 
inspired in recent years by the mode-coupling theory for idealized glass transitions 
(MCT).'2l Within this theory, the transition from a liquid to an idealized glass state 
is described by a bifurcation in the solutions of certain polynomials in a variable /: 
Upon smooth variations of a control parameter, say, the temperature T, variable / 
jumps from / = to a finite critical value f = > at T = Tc and increases further 
with the distance from the transition Tc — T. Variable / is defined by the long-time 
limit of some autocorrelation function (j){t); states with / > are identified as glass 
states, while states with / = shall be called fiuid or ergodic. 

The mentioned bifurcations originate from the equations of motion for density 
autocorrelation functions (j)q (t) of a system of N particles in a volume V with density 
p = N/V . The time-dependent fluctuations in the density are used to deflne the 
canonical normalized correlation function (j)q{t) = {p{t)* p) / {p* p) as is well-known 
in liquid-state theory.^' In this case, the correlation functions depend on the wave 
vector modulus q of the corresponding Fourier transform in space. Using projection- 
operator techniques, one can derive the following equations of motion 

dUq{t) + i^gdtMi) + ^gMt) + j ^t' mg{t - t')dt' (t>q{t') = , (Ma) 

with characteristic frequencies 17^, a white noise i/g, and the memory function 

1 I' d^k 

m,{t) = H<l^k{t),V] = - / , (Mb) 
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where the interaction potential is encoded in the vertex 

^q,k = P^'i^kS^j:P qkck + q{q-k)c^^_J^^ jq^ , (1-lc) 



2 ,4 



through the static structure factor Sq of the fluid and Cq its direct correlation func- 
tion. Both functions can be calculated from the interaction potential together with 
some closure relation that is known typically only in some approximation^ For 
t — 7- oo, one gets an algebraic equation for the correlators' long time limits (\)q{t) — > /g, 
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It was discovered in 1984 that Eq. (|l-ldp can exhibit nontrivial solutions, fq > 0, 
for microscopic interaction potentials and realistic values of the density in the sys- 
temP It was shown later that only singularities of type Ai can occur in Eq. (jl-ldp p* 
and these singularities are equivalent to those emerging for the parameter space of 
roots of polynomials upon variation of the coefficients.—' Hence, an A2 singularity - 
also called fold - signals a double root in the solutions like in the equation + 1 = 
for tc = 0. For a polynomial of high enough order, the variation of a single control 
parameter is sufficient to encounter an A2 singularity. Generically, the variation of 
i — 1 control parameters is necessary to identify singularities of type A^. An ^3 
singularity - called cusp - requires two control parameters to adjust a polynomial 
to a cubic root; an A^ singularity - the swallowtail - requires the variation of three 
parameters. It was shown in a theorem by Whitney that only A2 and A^ are robust 
singularities, all other singularities can be removed by small perturbations of the 
control parameters .'3' 'El 

Within MCT, the A2 singularity can be identified with a liquid-glass transition 
if / jumps from / = to a finite value at the transition. Once such a singularity 
is identified, asymptotic expansions of Eq. can be used to derive the long-time 
behavior of the correlator (j){t). For the A2 singularity, these asymptotic expansions 
yield two-step relaxation, time-temperature superposition, and power-law scaling.'^' 
The unique number characterizing the leading terms of the asymptotic expansion 
for an A2 singularity is the exponent parameter A which is between 0.5 and unity. 
In addition to liquid-glass transitions, fold singularities can also describe glass-glass 
transitions: In this case an existing first glass state with / > /f transforms into a 
second distinct glass state with / > /I > /f discontinuously. The endpoint of a 
line of glass-glass transition points is the ^3 singularity. An A4 singularity signals 
the emergence of a glass-glass transition line from an otherwise smooth surface of 
liquid-glass transitions. Every A^ singularity is characterized by a unique number 
IJ,e > that determines the properties of the asymptotic expansions.'^ /i£ — )• signals 
the emergence of the higher-order singularity A^^i where in turn /x^+i defines the 
leading terms of the expansion. An A2 singularity's exponent parameter is identical 
to /i2 = 1 — A. Therefore, a fold gives rise to a cusp once A approaches unity. 
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§2. Glass Transitions, Glass-Glass Transitions, and Higher-Order 

Singularities 

In the following, the transition singularities of MCT shall be reviewed briefly for 
hard spheres (sec. 12. ip . sticky hard spheres (sec. 12. 2p . and the square- well potential 
(sec. 12. Sp . The section [2] shall be concluded with the discussion of newly discovered 
transitions in the square-shoulder system (sec. 12. 4p where the static structure factors 
are evaluated in Roger- Young (RY) approximation. In sec. [3l new results shall be 
given for the square-shoulder system with the static structure factors calculated from 
the Percus-Yevick (PY) approximation. The comparison between results from RY- 
and PY-calculations gives insight into the robustness of the predicted glass-glass 
transition phenomena. 

2.1. Glass Transition for the Hard-Sphere System - a Fold 

MCT was first applied to the hard-sphere system (HSS) where a glass transition 
was identified upon varying the control parameter packing fraction ip = 7rpd^/6 
for particles of hard-sphere diameter d^ The transition is predicted for a packing 
fraction of ip"^ = 0.516 when using the PY approximation for the calculation of the 
static structure factor The hard-sphere interaction can be realized to a very 

good degree in experiments performed in colloidal suspensions. In such experiments 
the glass transition is found around a packing fraction ip^ = 0.58, moreover, the two- 
step relaxation, scaling laws and several other features of the theoretical predictions 
are confirmed.® From a thorough analysis of Eq. p-lcp one can derive that for 
the HSS, the glass transition is driven by the hard-core repulsion as encoded in 
the principal peak of the static structure factor. The HSS has been investigated by 
asymptotic expansionJ^ with the exponent parameter being around A = 0.7 which is 
close to the typical value for many other glass forming substances. The full evolution 
of glassy dynamics over eight orders of magnitude in time for various densities has 
been demonstrated for a tagged particle's mean-squared displacement (MSD).'^^ 

2.2. Glass-Glass Transition in the Sticky Hard-Sphere System - a Gusp 

For the investigation of an singularity, more than one control parameter needs 
to be varied. Generically the variation of two control parameters extends an isolated 
A2 singularity to a line of A2 singularities, so the existence of an A^ singularity is by 
no means certain. However, if such an A3 singularity exists in the given parameter 
plane, it can be found as an endpoint of a line of A2 singularities where close to 
the A^, these A2 singularities can be identified as glass-glass transition points. At 
the A^, the distinction between the different glass states, the discontinuity between 
the fq on both sides of the glass-glass transition, vanishes. Such endpoint singu- 
larities have been described first for so-called schematic models of MCT where the 
microscopic details, i.e., the g-dependence, has been dropped in favor of mathemat- 
ical simplicity. Rather than a precise interaction potential these models capture the 
mathematical structure of the problem. In these schematic models the transition 
points can be calculated analytically and the asymptotic expansions are not affected 
by limited numerical accuracy. For such schematic models, the asymptotic expan- 
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sions for singularities with i > 2 have been performed: Two-step relaxation and 
time-temperature superposition become invalid, and - most notably -~ logarithmic 
decay laws emerge P''^ 

The first microscopic model with an singularity in its parameter plane was 
found for Baxter's sticky hard-sphere model (SHSS).^^' The SHSS adds a short- 
ranged attraction to the hard-sphere potential similar to the square-well interaction, 
cf. Fig. [1] (a). However, within the SHSS, depth F and width 6 can only be 
changed together while the product r6 remains fixed, in addition the limit 5 ^ is 
performed; this defines a so-called stickiness parameter r. Hence, the SHSS has two 
control parameters, packing fraction ip like the HSS and the stickiness r. For the 
SHSS, the MCT predicts glass-glass transitions and an A^ endpoint singularity P^J'^^ 
The glass-glass transitions take place between a repulsion-dominated glass state - as 
known from the HSS - and an attraction-dominated glass state that resembles bond 
formation as known from the gelation transition.tHl' While for a treatment of the 
SHSS within MCT, a cutoff wave-vector space needs to be introduced, it is by now 
well-known how this cutoff can be interpreted as an inverse length scale, and how 
the MCT results stay well-defined.'^ In accordance with Whitney's theorem, for 
small changes in the cutoff neither fold nor cusp singularities change qualitatively. 
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Fig. 1. (a) Square-well potential with three control parameters: packing fraction (p = -KpcP /&, well 
depth r = uo/kBT, and well width 5 for particles of diameter d at density p. (b) Square-shoulder 
potential with shoulder height F — uo/ksT and shoulder width 5. 



2.3. Glass-Glass Transitions in the Square-Well System - a Swallowtail 

In order to avoid entirely the introduction of a cutoff and other peculiarities 
of the SHSS, one can extend the model attraction to finite widths in the square- 
well system (SWS), cf. Fig. [1] (a). Here, the control-parameter space becomes 
truly three dimensional, the parameter triple {ip, F, 5) defines each state. The MCT 
glass-transition scenarios have been worked out for the SWS,'^ and in addition to 
the cusp scenario of the SHSS there exists a characteristic well width 5* for the 
SWS above which no cusp singularity can be found in the {ip,F,6) parameter plane 
for fixed 6 when 6 > 6* . When 6 < 6* , the SWS always exhibits a glass-glass 
transition line with an A^ singularity as endpoint. For the exceptional point 5 = 6*, 
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the glass-glass transition lines together with the line of singularities vanish in 
an A4 singularity, giving rise to a three-dimensional geometric structure known as 
swallowtail. In contrast to the stable A2 and A3 singularities, the A4 is an isolated 
point in the parameter space that is sensitive to small numerical deviations from 
(99*, r*,6*). But while the details of the approximations involved and the numerical 
implementation change the location of the singularity, its existence is a robust 
prediction of MCT for systems with short-ranged attraction regardless of the specific 
interaction modePSi or closure relation for the static structure factors.^ 

The occurrence of the glass-glass transitions can be traced to two different mech- 
anisms of arrest - resulting from the interplay of repulsion and attraction - in the 
vertex, cf. Eq. (|l-lcp . The repulsion-dominated glass state within MCT is produced 
by the local structure on the wave- vector scale of the peak of Sq. The second mech- 
anism of arrest is given by a 1/g-tail in Sg for large wave vectors Once that 
tail has enough weight in the vertex, a second transition - the glass-glass transition 
- can occur. This possibility of an additional transition is robust regarding the clo- 
sure relation as long as the tail can become large enough. This independence of the 
results from the closure relation was demonstrated explicitly for the PY compared 
with the mean-spherical approximation .f^^'l^ 

The asymptotic solutions have been worked out in detail for the SWS and in- 
volve logarithmic decay laws,^ Vogel-Fulcher-like divergence of time-scales,'^ un- 
conventional (i.e., non-power-law) critical decays at the higher-order singularities ,123 
and the interplay of two A2 singularities at the crossing of glass- and gel-transition 
lines. Experimental verifications of the scenarios predicted for the SWS are found 
in computer simulations and in colloidal suspensions with attraction among the 
particles. Confirmations include the reentrant behavior of the lines of A2 singu- 
larities,'^"'^ the dynamics at a crossing of glass- and gel-transition lines for mi- 
cella^^ and colloidal suspensions,—* and the logarithmic decay of correlation 
functions together with novel power laws for the MSD."'^'^''^'^' 

2.4. Glass-Glass Transitions in the Square-Shoulder System 

When the attractive well of the SWS is replaced by a repulsive step, one obtains 
the square-shoulder system (SSS), cf. Fig. [1] (b). The SSS also has three control 
parameters, {ip, T, 5), with F now symbolizing the height of a shoulder. Allowing the 
control parameter F to carry a sign, SWS and SSS may even be plotted into a single 
diagram. Recently, the MCT predictions for the SSS have been worked out using 
the RY approximatioEp^ for SgW^ Further details and references to the numerical 
algorithms can be found in recent work .'23' 

Figure El shows the results of MCT for the SSS at 6 = 0.145: For F = 0, the 
glass transition line (lower panel) emerges from the HSS value of y^n^g = 0.5206, 
increases in density until around F = 3, bends over in an S-shape towards the HSS 
limit for the outer core, ip = 0.5206/(1-1-0.145)^ = 0.3468. The form of the transition 
curve can be understood in detail from the distortions of the local structures by the 
presence of the repulsion at distance 1-1-5 and the resulting changes to the principal 
peak of the static structure factor.'^S For both reentrant transitions it can be shown 
how the weakening of the local structure is compensated by either higher density (in 
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Fig. 2. Glass-transition diagram (lower panel) and exponent parameter A for the SSS using the RY 
approximation at 5 = 0.145. Diamonds show the transition points from an ergodic state to a 
glass state; filled circles exhibit the glass-glass transition points. Dotted lines show the limits 
for the hard-sphere system with diameter 1-1-5. 



the case of melting by cooling, i.e., melting at increasing F) or lower temperature 
which is equivalent to higher F (in the case of melting by compression). In the 
pressure-temperature representation, such a reentrant behavior can be identified 
with a so-called diffusion anomaly which is known experimentally e.g. for water.ISS) 
The exponent parameter A (upper panel) varies very little along the glass-transition 
line and stays around the common value of A = 0.7. From this unremarkable value 
for A no conclusion can be drawn about any glass-glass transitions or higher-order 
singularities in the vicinity. 

Moving further into the glass state, fq experiences an additional jump marked 
by the full circles in Fig. [2l In contrast to the situation for the SWS, for the SSS the 
glass-glass-transition line is located completely within the glassy state. This gives 
rise to two endpoints, two singularities. These two endpoints are seen also in 
the upper panel of Fig. [2] where A approaches unity on either side. This novel line 
originates from the competition of the two repulsive cores that causes a beating in 
the static structure factor for large wave vectors; this beating ~ if sufficiently large in 
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amplitude - can bring additional weight to the MCT vertex in Eq. ()l-lcp and hence 
cause an additional discontinuous transition.!^ In physical terms, the localization of 
the already arrested particles drops drastically when the glass-glass transition line is 
crossed since upon crossing the line, the localization now happens at the outer core 
rather than at the inner core as before. When varying the control parameters, the 
beating is most pronounced when inner and outer core have about the same impact 
on Sq - if either one of the cores is dominant, the beating is diminished and the 
discontinuous transition ends in two A3 singularities, respectively. Physically, the 
glass-glass transition becomes impossible at the upper endpoint when the density 
becomes too high for a transition to the outer core; the glass-glass transition also 
vanishes at the lower endpoint because the density becomes too low at the respective 
shoulder height to force the particles closer together. For shoulder widths 5 larger 
than the value shown in Fig. [2l the glass-glass transition line moves towards and 
merges with the glass-transition line.^^ 

§3. Square-Shoulder System using the Percus-Yevick Approximation 

In this section, it is demonstrated how the scenario shown in sec. 12.41 changes 
when the mentioned beating has not enough weight in the vertex (|l-lc|) to cause an 
additional transition line. 

Figure [3] shows the glass-transition scenarios for various shoulder widths of the 
SSS when the PY approximation is used for the calculation of the static structure 
factor. The overall behavior of the transition diagram is similar to the RY results: 
The glass-transition first increases in packing fraction ip when starting from the HSS 
limit (puss = 0.516 at = 0. Around F = 2.5 the curve bends downwards and 
reaches the limit of the HSS with an outer core of (/?hss/(1 + at around F = 5. 
Different from the RY result, the S-shape of the transition curve for intermediate 
values of F only develops at higher values for the width 5. Also in contrast to the 
RY result, there appears to be no indication of additional discontinuities in the fg 
or endpoint singularities inside the glass regime. However, while A in the upper 
panel of Fig. [2] stays almost constant at the glass-transition line for the RY results, 
for the glass-transition line within PY approximation the exponent parameter A 
varies considerably with F for any given width 5 in the upper panel of Fig. [3l In 
addition, the shape of the A-versus-^T curves varies drastically with 6. For 5 = 0.25, 
an exponent-parameter maximum is around A = 0.85 which is already rather high. 
For larger 5, the A-versus-^T curves exhibit double maxima over relatively small 
intervals in F, like between 1.85 and 2.1 for 5 around 0.28. For 5 = 0.2785 and 
6 = 0.28 the parameter A approaches unity very closely at the right and the left 
end of the interval, respectively. Beyond that regime, A decreases again for larger 6, 
and as shown for 6 = 0.3, the separation of both A maxima increases; the A maxima 
are located at -T = 1.65 and F = 2.4, respectively, while their values are still as 
high as 0.9. In summary, while no glass-glass line can be detected, the exponent 
parameter approaches A = 1 very closely. These results indicate a nearby higher- 
order singularity without the presence of glass-glass transitions. More puzzling still 
is the existence of two maxima in A very close to each other. 
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Fig. 3. Glass-transition scenario (transition lines in the lower and exponent parameter A in the 
upper one) for the SSS using the PY approximation for 5 — 0.25 (dotted line), 0.2785 (full line), 
0.28 (dashed line), and 0.3 (chain line). In the lower panel, the curves for S = 0.2785 and 0.28 
are almost on top of each other. 

To clarify the nature of the higher-order singularities on the glass-transition 
curves, the glassy region is inspected in more detail for 6 = 0.25 in Fig. [H The 
behavior of A indicates that here the higher-order singularities are further away from 
the transition line, but a rather broad maximum already hints at these. It is well- 
known that after crossing an A2 glass-transition singularity, the long-time limits 
fq increase above their critical value in a square-root in the control parameter; 

i.e. in the HSS, this increase is proportional to y^ip — '/^hss-'^ When a glass-glass 
transition occurs, this square-root increase after the first transition is superseded 
by an additional discontinuity in fg followed by the square-root increase after this 
second transition. While this additional discontinuity was used to identify the glass- 
glass transition line in Fig. [51 such discontinuity is absent for the PY calculations. 
Nevertheless, the evolution of the fq within PY approximation shows a characteristic 
deviation from the square-root behavior in certain parameter regions: First, the 
range of validity of the square-root increase is sometimes far smaller than known from 
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the HSS; second, the apparent square-root increase at larger distances from the glass- 
transition indicates a square-root with a different extrapolated transition point c^'^pp 
than the one given by the discontinuity in fq at (p'^. This difference between the actual 
and the apparent transition point can be quantified as a difference in the control 
parameters, say as a relative difference in packing fraction, A(p = (c^^pp — (p^)/Lp'^. 
When this relative difference in packing fraction exceeds 1%, this anomaly is marked 
by open circles in Fig. 21 It is seen that this line of fq anomalies strongly resembles 
the line of glass-glass transitions in Fig. [2l It can therefore be concluded that within 
MCT both approximations, PY and RY, yield similar glass-transition scenarios with 
a possible line of glass-glass transitions inside the glassy regime. 

§4. Conclusion 

In the present work, the glass-transition diagram for the SSS has been calculated 
for the PY approximation. These results can now be used to estimate the robust- 
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ness of the results for the SSS obtained with a different closure relation, e.g. with 
results from the RY approximation.'^S While the comparison of different closure 
relations for the SWS only resulted in shifts of the control parameters, the situation 
is more involved in the case of the SSS. The existence of the disconnected glass-glass 
transition line depends on two trends that both become more prominent for larger 
values of 6. First, the vertex in Eq. (jl-lcp obtains additional weight for higher wave 
vectors through a beating phenomenon as described above; this makes a glass-glass 
transition possible if the weight becomes strong enough. Second, at the same time 
this supposed glass-glass transition line moves closer to the glass transition line and 
merges with it. For the RY approximation, the glass-glass transition line becomes 
manifest through a discontinuity in fq and then moves towards the glass transition 
line. In contrast for the PY approximation, the not-yet-manifest glass-glass tran- 
sition moves towards the glass transition line and merges with it before it is fully 
developed as a discontinuity in the fq. When the endpoints of this hidden line of 
glass-glass transitions merge with the regular glass transition line, very high values 
of A result since the first trend of increased weight for a glass-glass transition keeps 
increasing. In this sense, the glass-glass transition line emerges extremely close and 
on top of the glass transition line and is at the same moment absorbed by the glass 
transition line. 

The differences between RY and PY approximation with respect to the higher- 
order singularities and the corresponding glass-glass transitions are highly non-trivial 
in the theoretical calculations. However, for the experimental test if such a scenario 
exists, both scenarios need to be observed in combination: Numerical deviations 
similar to the ones shifting the HSS-PY transition from ip'^ = 0.516 to 0.58 in the 
experiment can for the SSS switch from the RY scenario shown earlier*^ to the PY 
scenario described here. 
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